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Abstract 


Hadron therapy is a novel treatment against cancer. The main advantage of this therapy 
causes less side effect in comparison to X-ray irradiation methods. Hadron therapy is just 
ahead of a significant breakthrough since this technique can be more precise, applying proton 
computer tomograph (pCT) to map the stopping power in the tissues. 

The research and development of a pCT require a fast detector to measure the energy of 
hadrons behind the patient. The best detector option is called hadron-tracking calorimeter, 
which consists of sandwich layers of silicon tracking detectors and absorber layers. The combi¬ 
nation of measuring the trajectory (tracking process), and, in parallel, the energy of relativistic 
particles, can provide high-resolution hadron imaging. This semiconductor-based technology 
requires stable temperature and homogeneous cooling. 

1 have worked in the development of this detector in the Bergen pCT Collaboration for two 
years. Last year my work was to investigate the temperature distribution in the calorimeter and 
examine two cooling concepts in detail. 1 performed both analytical and numerical calculations 
to analyze the temperature distribution of the calorimeter. The final decision about the design 
takes into account many engineering aspects, such as reliability, flexibility, and performance. 
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1 Introduction 


1.1 Cancer 

Cancer has become one of the leading reasons for death in the developed world. It is 
responsible for 25% of all death in Hungary [1], It affects all age groups, but the risk of cancer 
is increasing with age, as one can see in Figure 3. 



Figure 1: Cancer affects all age groups, but the risk of cancer is increasing with age. This 
statistic was made in Scotland in 2016 [2], 

This section is about the nature of cancer, briefly summarized based on the literature [3, 4, 5]. 

The cell cycle is an important period of life of the cells because they divide and replace 
dead cells. In healthy tissues, there are stem cells, which can divide at any time. However, for 
most of the cells, the number of divisions and their times are limited. This is not true for the 
cancer cells because a DNA mutation turns off this limitation. 

Usually, a cell dies after a lifetime, or in case of serious injury in its DNA. This function is 
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called apoptosis, which, in general, does not work in cancer cells. 

There are two types of cancer. The first is called bengilus. In the case of bengilus cancer, 
the tumor does not reduce the life expectancy and quality of life of the patient. The second 
one is called malignus, which means risk in life or the life quality of the patient. It is hard to 
determinate whether a tumor is bengilus or malignus. Usually, bengilus tumors expand slower, 
and they have a tegument, but the lack of these effects does not evidence that the tumor is 
malignus. 

Unfortunately, tumors can expand in the surrounding tissues. This phenomenon is called 
an invasion. Cells of the tumor can spread far from the original one, this is called metastasis. 
After the treatment of cancer, it is still possible to find other cancer cells in the surrounding 
tissues that can cause a new, so-called residual tumor. 

There are status groups, called stages, to categorize the tumors. The traditional staging 
has five stages: 

0. stage: there is no invasion, 

I. stage: small, local tumor, with small invasion (without metastasis), 

II. stage: local tumor with more significant invasion, probably with metastasis in the local 
lymph nodes, 

III. stage: local tumor with significant invasion with high probability to has big metastasis in 
the local lymph nodes, 

IV. stage: local tumor with huge invasion and with far-reaching metastasis. 

1.2 Mapping methodes in the detection and treatment of cancer 

There are lots of methods to measure or find a tumor in the human body. This subsection 
presents them briefly based on [6]. 

Ultrasound method uses sound waves to take a picture from the inside of the body. Applying 
this method, the sound wave penetrates into the body of the patient and is scattered on the 
boundaries of the organs with different speed of sound. A microphone records that reflection 
and a computer reconstructs the map. A usual ultrasound machine takes a 2D picture from 
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a segment of the human body, but nowadays it is also possible to obtain a 3D picture. The 
ultrasound technology does not have any side effect. 

An X-ray computer tomograph uses X-ray beams to fluoroscope through the human body. 
It creates a 3D picture from the nucleon density of the body. This type of imaging causes 
radiation in the body, which can damage healthy tissues. 

The MR technology is based on the measurement of the magnetic held caused by the water 
molecules in the human body. The water molecules are excited by a strong magnetic held 
generated by superconductor magnets. The detector measures the changes of the magnetic 
held after the magnets are turned oh. For instance, it makes visible the blood how in the veins. 
MR technology does not cause any side effect, based on the current knowledge of science. The 
mapping process takes a long time, and the patient have to lie in a narrow, and closed space 
that causes feeling of being locked up in 4 - 6% of the patients. 

Mapping with radiative isotopes can give information about the life cycle of a tumor. This 
technology is based on the density of the chemical material which is absorbed by the tissues 
and the tumor. Usually, the tumor absorbes different amount of chemical element than the 
surrounding healty tissues. Thus, the concentration difference of the isotopes is measured here. 
The directional dependence of the radiation makes possible to construct the 3D map. Moreover, 
the amount of the absorbed radiation offers information about the type, size and the current 
state of the tumor. It is also useful to find metastasis; however, the irradiation of healthy 
tissues is unavoidable. 

1.3 Treatment of cancer 

The treatment aims to ensure high life expectancy (or as high as possible) and life quality (or 
as high as possible). In order to reach this goal, modern health care uses surgery, radiotherapy, 
chemotherapy, and hormone therapy. These treatments can be applied either standalone or 
combined as a complex therapy to improve the effectiveness. That subsection is summarized 
based on [7, 8, 9]. 

Every treatment starts with canceling the tumor. It is crucial not to leave any cancer cell in 
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the surrounding tissues because it can cause a residual tumor later. If it is not possible due to 
some reason, then it becomes necessary to kill these cells in the second part of the treatment. 
That part is about to kill all the remaining cancer cells in order to avoid metastasis. These 
cells are probably to be a metastasis, or cells in the vascular and lymphatic system, which try 
to find a new living space and create a new metastasis. 

Surgery is the leading treatment of cancer. It is the oldest treatment of tumor diseases. In 
general, it is used to remove the original tumor, but it is useful in removing the metastasis, too. 
Its one of the most significant advantage is that the cancer cells cannot become resistive to the 
surgery. During this process, doctors cut out the tumor and usually remove the local lymph 
nodes also from the body to reduce the risk of metastasis. The surgery burdens the body of the 
patient. If the patient is not in adequate conditions, the operation becomes impossible. There 
are other situations in which the outcome is the same: either the tumor is in a hardly available 
location, or too big to cut out. In these cases, doctors have to find another type of treatment 
or use complex methods. 

Radiotherapy uses the effects of the ionization of X-ray or ion beams. The beam exerts its 
effect in two ways, directly and indirectly. The direct way is when the particles of the beam 
damage the DNA and the critical parts of the cell. The indirect way is when the particles of the 
beam ionize the water molecules, and these molecules damage the DNA of the cancer cells. In 
both cases, the damage induces a self-destruct function of the cancer cell, thus the cell destroys 
itself. This self-destruction process is called apoptosis. Unfortunately, there are such cancer 
cells which have resistance to the radiation, the apoptosis process is prevented. 

Usually, chemotherapy uses apoptosis to destroy cancer cells. In the case of chemotherapy, 
doctors give medicine to the patient, which consists of chemical elements, inducing the apoptosis 
process. The most common chemicals are called cytostatics. The cytostatics affect the cancer 
cells more seriously than healthy cells. Besides, there are other chemicals to reduce the effect 
of cytostatics in the healthy regions or increase their exterminating effect in the cancer cells. 

Surgery as a standalone treatment is usual for a tumor in the I. or II. stages. In the III. stage, 
radiotherapy is recommended. Chemotherapy is typical for IV. stage. 
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Beyond these possibilities, sometimes the surgeon must remove the entire organ, this is 
called ‘organ removal surgery’. However, this is not preferred and tried to be avoided. For 
instance, using radiotherapy first reduces the size of the tumor, hence it is possible to preserve 
the organ as much as possible. Moreover, chemotherapy helps to reduce the risk of metastasis 
with canceling the remaining cells in the vascular and lymphatic system. 

1.4 Radiotherapy 

Radiotherapy uses the ionizing effect of radiation to destroy cancer cells. There are different 
ways to transfer the radiation to the tumor [10]. As a first option, one can use beams, which 
penetrates the tissues and the tumor, and in the meantime, it damages them. This treatment 
is called external beam radiation therapy. Its advantage is that it is usable in many segments of 
the body, without using surgery. Secondly, one can situate a radioactive source in the vicinity 
of the tumor, for example, into a body cavity. This technology is called brachytherapy. The 
advantage is the same, no operation is needed. As a third option, it is possible to ingest certain 
radioactive isotopes into the body. In fact, a tumor absorbs more isotopes than any of the 
other healthy organs. This treatment is called systemic radiation therapy. The disadvantage 
of this technology is that the patient becomes radioactive for a short time, so he or she has to 
stay in the hospital for some clays. 

One of the primary development in radiotherapy is to reduce the side effects caused by the 
radiation in healthy tissues. If one improves the ratio of the ionization in the tumor, it reduces 
the amount of ionization in the healthy tissues, thus it reduces the side effects. 

X-ray photon beams are the most common in external beam radiation therapy. The usage 
of protons or heavier ions instead of X-ray photons results in less ionization in the healthy 
tissues [10, 11]. Its advantage is based on the radiation concentration distribution. The X-ray 
photons generate the highest radiation after the entry point in the body, then the radiation 
is monotonously decreasing. Using ions instead of X-ray photons, it is possible to concentrate 
the radiation onto the tumor, since the ions stop in the tumor mostly. Ions create the highest 
radiation right before they stop. One can see the resulted radiation by X-ray photons and 
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protons in Figure 2. One can also see the radiation of each treatment in a cross-section of a 
head in Figure 3. 



Figure 2: The energy-loss of X-ray photons and protons in function 
of the amount of material in their previous path [12]. 
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Figure 3: Radiation in a head in case of proton and X-ray therapy 
[ 12 ]- 

As it is apparent, the application of ions in the external beam radiation therapy is beneficial. 
Nowadays, the number of hadron therapy centers is increasing. One can see the treatment room 
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of a hadron therapy facility in Figure 4. However, there is a design barrier in the spread of 
hadron therapy centers. This barrier is the lack of an accurate imaging tool for dose planning. 
Doctors use X-ray computer tomographs (CTs) to take a three-dimensional map from the 
nucleon density of the body. However, the energy loss of ions depends on the electron density 
of the body. It is possible to calculate the electron density using the measured nucleon density, 
but there is no simple relation between them, because of the different ratio between proton 
and nucleon number of atoms of the human body. As a result, one can calculate the electron 
density with 1.7% statistical error from nucleon density [13]. If one could measure the electron 
density map of the human body directly with the usage of ions, the statistical error could be 
reduced to 0.5% without significant error [13]. Overall, the development of a hadron computer 
tomography (usually called proton computer tomography or pCT) is beneficial. 



Figure 4: The patient is fixed into a bed during the radiotherapy. 
The patient is moved with the bed by a robot arm that can be 
moved in all directions. First of all they take a 3D map of the 
tissues nearby the tumor of the patient, after they treat the tumor 
by a proton beam [14]. 
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1.5 Proton computer tomography 


The original idea of computer tomography (CT) comes from Allan M. Cormack (1963) [15]. 
He won the Nobel Prize with this idea in 1979. 

A CT equipment moves around the body and takes pictures from the inside in several 
directions, and a computer calculates the 3D map using these pictures. In case of proton CT, 
one can use a beam to fluoroscope through the body. One scans the body with this beam in 
a particular direction, so one obtains a two-dimensional picture about the body. After, the 
detector is turned into a different direction, and the scanning process starts again. One can 
see the scanning and rotating process in Figure 5. Finally, a computer algorithm calculates the 
three-dimensional map of the body from the two-dimensional pictures. 

Direction of rotation 
of the detector and the beam 



Direction of beam Patient 



Figure 5: Left: it is visible how the proton beam scans the body. 
Right: it is visible how the detector and the beam rotate around 
the patient. 


During the scanning process, one has to use another detector at the other side of the 
patient in order to measure the energy and the path of particles coming from the beam. Such 
a detector is called a tracking calorimeter since it measures both the energy and the tracks 
of hadrons. 10 7 - 10 9 ^ particle rate is necessary to obtain a three-dimensional picture 
within a reasonable time for clinical use [16]. To measure the path of the particles with such 
a rate was impossible five years ago. However, thanks to the development of tracking silicon 
pixel detectors in the last years in CERN LHC, there are available tracking detectors that can 
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meet with this particle rate now [17]. These tracking detectors can measure the path of the 
particles, but there is not any available detector that can measure the energy of each particle 
with this particle rate. 

The Bergen pCT Collaboration aims to develop a calorimeter (calorimeter is a detector, 
which measures the energy of hadrons) build of alternating tracking detector and aluminum 
absorber layers, as a sandwich structure [16]. One can see the concept of this detector in Figure 
6. As a tracking calorimeter, we are going to use a tracking detector, which was developed in 
ALICE, CERN LHC, and which is called ALPIDE. 



My role in the Bergen pCT Collaboration is the investigation of the temperature distribution 
in the calorimeter and compare two cooling system concept. That work is crucial to ensure the 
accuracy of the detector. That part of the detector development is clearly separated from the 
other tasks related to the design. Thus the present thesis reflects my work. 

1.6 ALPIDE tracking detector 

ALPIDE was developed to upgrade the Inner Tracking System (ITS) of ‘A Large Ion Collider 
Experiment’ (ALICE) in 2019-2020, which is one of the four large experiment programs of the 
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Large Hadron Collider of the European Organization for Nuclear Research (CERN) [18]. One 
can see the ALPIDE detector in Figure 7. 



Figure 7: ALPIDE detector [19]. 


ALPIDE is a monolithic active pixel sensor (MAPS) type silicon detector [17], which con¬ 
tains the sensitive sensors together with the readout electronics on the same silicon layer, see 
Figure 8 for details. When a charged particle goes through the sensitive (light blue area, called 
Epitaxial Layer P-) area of the detector, it generates electron-hole pairs. The electrons and 
the holes are transferred in the sensitive area by diffusion. The holes are absorbed by the 
Substrate P++ layer (dark blue). There is an area (white), which is evacuated by electric 
voltage. When an electron reaches this zone, the electric field of the evacuated zone transfers 
it into the NWELL diode (green). This diode collects the electrons, then a microelectronic 
amplifier amplifies the signal. Another microelectronics decides whether the signal is higher 
than a threshold or not. If so, the microerectronic device sends a signal to the readout electron¬ 
ics with the coordinates of the pixel. This is a state of the art solution because the sensitive 
silicon layer contains the complete amplifier electronics. It is possible, thanks to the DEEP 
PWELL (red), which separates the NWELL parts of the microelectronics, thus the NWELL 
parts cannot operate as an undesirable diode. 
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Figure 8: Main parts of a monolithic active pixel sensor 
in a section view of one pixel [17]. 


The electrons, generated by the same particle, can be transfered to more than one pixels. 
This effect causes that more than one pixel sends signals to the readout electronics. We call 
as cluster the pixels send signals after a particle go through the detector. The size of the 
cluster depends on the absorbed energy in the detector, so it is possible to use this information 
to increase the accuracy of the calorimeter. However the cluster size also depends on the 
temperature of the detector, so it is very important to know the accurate temperature of the 
detector, because a small difference can cause error in the energy measurement. 


13 




















2 Detector 


2.1 Concept of the detector 

Let us recall the concept of the detector in Figure 9. The detector is built of two parts. The 
first part is a tracking detector, made of two layers of ALPIDE sensors. This part measures 
the direction of incoming particles. The second part of the detector is a calorimeter. It has 
a sandwich structure, and it is made of 35 ALPIDE sensor layers, separated by 4 mm thick 
aluminum absorber layers. This calorimeter measures the energy of the particles. 



One can see the structure of a detector layer in Figure 10. The ALPIDEs are glued on 
both sides of the aluminum absorber layer, alternately with a small overlap, hence the absorber 
functions as a support structure also. Here, two definitions are introduced for the layers, 
depicted in Fig. 10. Its first meaning is about the engineering aspects, and consists of the 
aluminium absorber with APLIDE sensors on its both sides. Its second meaning is about the 
data analysis aspects, which considers the ALPIDE sensors between two aluminium absorbers. 
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Figure 10: This figure shows the structure of a layer. Gray: ab¬ 
sorber, red: ALPIDE, blue: chip cable and orange: glue. 



Figure 11: A stack is a group of ALPIDEs, which is connected to 
the readout electronics with the same chip cable. 
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The ALPIDEs are grouped into stacks that form a readout unit. One can see a stack on 
one side of the absorber in Figure 11. The stacks are placed in this way, but they are mirrored. 
The thickness of the glue, ALPIDE, and chip cable are small (Figure 10) but not necessarily 
negligible. The resulted temperature difference AT in the perpendicular direction of the layer 
can be calculated with the following equation in steady-state [20]: 

AT=^, (1) 

where q is the heat flux, A stands for the thermal conductivity, and h being the thickness. 
One can find the constants and the results in Table 1. Indeed, these AT values are small, 
thus in later calculations their thermal resistances are neglected. A layer model is depicted in 
Figure 12. The volumetric heat generation occuring in the ALPIDE detectors is modelled as 
an equivalent volumetric heat generation in the absorber layer. 


Layer 

Material 

Thermal 

Thickness 

Typical 

Temperature 



conductivity 


heat flux 

difference 



A 

r w i 

mK 


h [/im] 

q 

rwi 

_rrP. 


AT [K] 

Glue 

Glue 


0.22 

10 


410 

0.0186 

ALPIDE 

Silicone 


149 

50 


410 

1.38 x 10“ 4 

Chip cable 

Polyamide 


0.24 

50 


410 

0.0854 


Table 1: Thermal resistance parameters of ALPIDE, chip cable and glue. 

One can see the essential geometrical parameters and the thermal load of a real layer in 
Table 2. It is not worth to model it in its original arrangement. Consequently, a simplified 
model is prepared using the parameters in Table 3. There are several layers side by side in 
which almost equal heat generation occurs. Thus, applying a periodic boundary condition, one 
can define a ‘unit’ in which half of the air gap appears on both sides together with one layer. 
Between the layer and the air heat convection occurs. 
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Figure 12: Comparing the real detector to the model. Q v is the real 
volumetric heat generation in ALPIDEs, and Q V2 is the equivalent 
heat generation in absorber layer. 


Part 

Material 

Height 

Width 

Thickness 

Heat generation 



[mm] 

[mm] 

[mm] 


[W] 

rrrh 


Absorber 

Pure aluminum 

200 

300 

4 




ALPIDE 

Silicone 

15 

30 

0.05 

8.2 

10 6 

Chip cable 

Polyamide 

27 

285 

0.05 




Glue 

Glue 

27 

270 

0.01 




Air gap 

Air 

200 

300 

1 





Table 2: The essential parameters of the real layer. 


Part 

Material 

Height 

Width 

Thickness 

Heat generation 



[mm] 

[mm] 

[mm] 


W] 

.m3. 


Absorber 

Air gap 

Pure aluminum 

Air 

200 

200 

300 

300 

4 

1 

8.3 • 10 4 


Table 3: Overall dimensions and the heat load in the layer model. 
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2.2 Without any active cooling system 


Here, I consider the detector that is covered by a casing, it protects the detector from 
external effects. As a first approximation, I neglect the inner temperature difference, and 
the entire detector is modeled as a lumped capacitance. Despite this simplification, it offers 
quantitative information about the characteristic heating of the device, without any active 
cooling system. One can see the material data in Table 4, which are used to characterize 
the detector. In Table 5, further parameters are summarized for modeling. Table 6 contains 
the number of ALPIDEs and the heat generation is specified for each part. The amount of 
heat generation is determined based on the description of ALPIDE elements, provided by the 
research group. One can see the heat generated by the particles of the beam in Table 7. The 
heat generated by the beam is several orders of magnitude smaller than the heat generated by 
the electronics, so I neglected it. Regarding the heat transfer coefficient, I assumed air at rest 
around the detector casing. 


Part 

Density 

Mass 

Specific heat 

Heat capacity 



'kg' 

_mT 


[kg] 


r j 

[kg-Kj 



' J' 
.K. 


One layer 

2850 

0.684 


434 


297 

35 layers 

2850 

23.94 


434 

10395 


Table 4: Thermal parameters for the lumped capacitance model. 


Height 

Width 

Thickness 

Surface area 

Heat transfer coefficient 

[mm] 

[mm] 

[mm] 

[m 2 ] 


r w ] 

_m 2 -K_ 


200 

300 

175 

0.295 

5 


Table 5: The dimensions of the detectors, and the applied heat transfer coefficient. 
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ALPIDE 

Stack 

Layer 

Detector 

Number of ALPIDEs 

1 

9 

108 

3780 

Generated power 

185 mW 

1.66 W 

19.9 W 

697 W 


Table 6: The generated heat and the number of ALPIDEs. 


Energy of a particles 

Particle rate 

Generated power 

[MeV] 

j- particle j 

[mW] 

200 

10 9 

35.2 


Table 7: The generated heat by the particles of the beam. 


The lumped capacitance method is the mathematical formulation of the I. law of thermody¬ 
namics, i.e., one has to formulate the balance of internal energy. On the left hand side, there is 
the time evolution of the internal energy, considering constant specific heat c and mass density 
p. On the right hand side, I consider the thermal load and the heat transfer by convection [20]: 

dT 

— • m • c = Q- (T-Tqo) • a ■ A . (2) 

The solution of the differential equation (2) reads: 

T(t) = C cons t • e m-c* + Tqo H tA (3) 

a ■ A 

in which I can define the time constant: 


r = 


m • c 
a ■ A 




and the C cons t is about to consider the initial condition 


( 4 ) 


T(0) = T 0 . 


One obtains the following form of the T(t) function (3): 

T(t) = -^ + T oo + (T 0 -T oo -- ( ^-)-eV . 
a ■ A a ■ A 


( 5 ) 

( 6 ) 
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The maximum allowable temperature for the detector is 40 °C, which can be exceeded easily 
in continous operation, even for a one hour interval (see Figure 13), and Table 8 summarizes 
the constants used in the calculation. The allowable 40 °C exceeded under 230 s, which time 
interval is too short for the detection. Moreover, one should wait too much time to use the 
detector again. Consequently, an active cooling system is desired. 



Figure 13: The temperature history of the detector without active 
cooling. 


To 

Too 

Q 

a ■ A 

m ■ c 

r 

[°C] 

[°C] 

w 

rwi 

K 

f kJ 1 

K 

[s] 

25 

25 

697 

1.48 

10.4 

7030 


Table 8: The constants used to determine the temperature history. 
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2.3 Cooling concepts 


Here, I compare three significantly different cooling concepts: 

• water cooling at the edge of aluminum layers, 

• air cooling between the layers, 


• water cooling inside the aluminum layers. 


Two of them use water as coolant and one concept uses air. Let me review them in the following. 

A) Water cooling at the edge of aluminum layers. In concept A, the detector layers 
are cooled at their top and bottom edges. The scheme of this concept is depicted in Figure 14. 
The benefit, compared with concept C, is that it does not need any extra material between the 
layers, thus it does not cause any extra scattering, consequently, it does not have any negative 
effect in the data reconstruction. 


Water cooling 
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Water cooling 


Water cooling 
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"O 
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A-> 

m 

o 

to 


Water cooling 


Figure 14: Schematic arrangement of concept A. 


On one hand, this is a compact and quiet arrangement, so it is perfect to use in treatment as 
well. It does not need a high coolant volume flow because of the high density and heat capacity 
of water. One needs a heat pump to cool the coolant that can be set far from the detector, 
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and it makes possible to change the mean temperature of the detector. On the other hand, 
we cannot influence the temperature distribution. That could be a problem when a significant 
temperature difference occurs in a layer. In this case, one cannot use this concept of detector 
cooling. It reveals that the knowledge about the temperature distribution could be essential. 

B) Air cooling between layers. In concept B, there is forced airflow between the detector 
layers, which cools the detectors directly. One can see this cooling concept in Figure 15. Using 
this cooling method, almost the entire surface of the detector is exposed to forced convection. 
Here, I expect a more advantageous temperature distribution with less significant temperature 
gradient within a layer. 



Figure 15: Schematic arrangement of concept B. 

It could be easier to realize this concept since one can use a simple fan instead of a water 
circuit system, thus the engineering design is easier than the others. Moreover, it becomes 
completely irrelevant about how to protect the electronics from the coolant. However, it uses 
the air of the room as a coolant which restricts the minimum temperature and the cooling rate. 
It may result in high airspeed and disadvantageous vibrations in the structure. It is also an 
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interesting question to investigate the temperature distribution in this case. 

C) Water cooling inside aluminum layers. This concept can produce the most ho¬ 
mogenous temperature distribution. However, it is disadvantageous for the data analysis due 
to the flowing medium among the layer. Moreover, it would be significantly harder to protect 
the electronics in case of any failure. Later, this concept is excluded from further investigations. 
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3 Comparing the concepts 


Now, my aim is to present the basic differences between concepts A and B by comparing 
the resulting temperature distributions to each other. In the following, I consider only one¬ 
dimensional heat conduction along the layer (direction of x axis). In case of concept A I also 
used symmetry boundary condition at the middle (x = 0) of the layer. 

3.1 Concept A 

In general, the heat transfer coefficient can be high enough to replace the convection bound¬ 
ary condition with a so-called first-type boundary in which the temperature is prescribed di¬ 
rectly. This is what we apply here. Integrating the Fourier heat equation in Cartesian coordi¬ 
nate system, one obtains [20] 

t w-| x 2 + |(y) 2+t »- < 7) 

in which Tq is the coolant temperature that used as a boundary condition, and Li stands for 
the height of the detector. Due to the symmetrical arrangement, the maximum temperature 
occurs at the middle, i.e., considering x = 0, 

Tmax = lx (y) +t °- (8) 

The desired operation condition requires as small as possible temperature difference inside a 
layer, thus the difference between T max and Tq has the importance: 

T di ff = T max - T 0 = f A) . (9) 

One can see the temperature distribution in Figure 16, and the parameters can be found in 
Table 9. 

In the current state of the detector design, the research group did not decide the material 
of the absorber. Most probably, it will be made from an aluminium alloy. Previously, the 
material properties of pure aluminium is used. However, I have found important to investigate 
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Figure 16: Temperature distribution in the layer, where the direction of x axis is visible in 
Figure 14. 


Constant 

Notation 

Value 

Coolant temperature 

Thermal conductivity of the absorber 

Height of the detector 

Equivalent volumetric heat generation in the absorber 

To 

A 

Li 

q v 

25 °C 

237 W 
Z,il m-K 

0.2 m 

8.3 • 10 4 

in' 3 


Table 9: Constants used to the calculation of the temperature distribution in concept A. 

the possible temperature difference Tjjff when the thermal conductivity is changed, accounting 
the fact that the material might be different in the real design. Thus thermal conductivity 
from 100^- to 237^^ is considered, which are the realistic values for an aluminum alloy. The 
maximum temperature has a strong dependence on this material property as one can see in 
Figure 17. Most likely, the material will be A11050, which has a thermal conductivity of 222 ^^ . 
The difference in maximum temperature between pure aluminum and A11050 is 0.12 °C, which 
is small, and the calculations with pure aluminum are also a good approximation for A11050 in 
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steady-state situations. 



Figure 17: Maximum layer temperature as a function of thermal conductivity of the absorber 
layer. The temperature of the edge of the detector is set to 25 °C. 

3.2 Concept B 

3.2.1 Heat transfer coefficient 

The calculation of the heat transfer coefficient is not evident in this case. There is a wide 
and long, but very narrow gap between the layers. Considering it as a tube with an equivalent 
tube diameter does not work evidently. The heat transfer coefficient in narrow annular gaps 
was measured previously by Tachibana and Fukui [21]. Now, I am considering the equivalent 
dimension between a narrow annular gap and a gap with rectangular cross-section by wrapping 
the rectangular one around a cylinder (see Figure 18). One can find the geometrical results 
in Table 10. The equivalent tube diameter is necessary for the calculation of Reynolds (Re) 
number. It is calculated based on the dimensions of the annular gap, but one would obtain 
almost the same results using the rectangular cross-section directly. 

The general strategy to determines the heat transfer coefficient requires an empirical (or 
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Figure 18: Transformation between a rectangular cross-section gap 
and a narrow annular gap. 


Dimension 

Notation 

Calculation 

Value [mm] 

Hight of the air gap 

Li 


- 

200 

Width of the air gap 

l 2 


- 

300 

Thickness of the air gap 

l 3 


- 

1 

Equivalent outer diameter 

Di 


¥ + L 3 

96,5 

Equivalent inner diameter 

d 2 


¥-l 3 

94,5 

Equivalent tube diameter 

D e 

(D?-D2)tt 

(Di+D 2 )tt 

= Di - D 2 = 2 L 3 

2 


Table 10: Geometry of the air gap between two layers and the equivalent annular air gap. 

semi-empirical) relation between the Nussclt (Nu) number (the dimensionless heat transfer 
coefficient) and some other dimensionless numbers such as the Reynolds (Re), Prandtl (Pr), 
Grashof (Gr), depending on the particular situation. 

Here, 1 use the following relation to determine the Nusselt number [21]: 

Nu = 0.017 ^1 + 2.3^j ^ Re°' 8 Pr 1 / 3 , (10) 

where the Reynolds number and the Prandtl number are calculated using the equations [20]: 
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( 11 ) 


and 


v 


( 12 ) 


Re = 


wD f 


Pr = - 
a 


Using the Nn number, I can determine the heat transfer coefficient as 

NuA 


a = 


Dr 


(13) 


For this particular situation, one can see the heat transfer coefficient as a function of airspeed 
in Figure 19, and the constants are summarized in Table 11. 



Figure 19: The heat transfer coefficient as a function of airspeed. 


3.2.2 Temperature distribution 

In order to obtain the temperature distribution, I have applied the following simple numer¬ 
ical approximation and discretization. The solution is based on the internal energy balance, 
accounting the heat transfer among each cell. The aluminium plate and the air gap are splitted 
into n parts, see Figure 20 for details. Every cell of the aluminum (‘w-type’ cell) are connected 
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Constant 

Notation 

Value 

Unit 

Air temperature 

Density 

Kinematic viscosity 

Specific heat 

Thermal conductivity 

Thermal diffusivity 

T 

P 

V 

c 

A 

a 

25 

1.19 

1.58 • 10“ 5 

1013 

0.0261 

2.17- 10” 5 

u 

[°C 

'kg' 

_mT 

m 2 

s 

J 

cg-Kj 

W 

m-K 

m 2 

s 



Table 11: Parameters of the air. 


by heat conduction and connected by heat convection to the air (‘a-type’ cell). Heat con¬ 
duction among the a-type cells is not considered. It is possible to determine the steady-state 
temperature distribution by iteration among the cells. 


wO 

wl 

w2 

... 

wn 

aO 

al 

a3 

... 

an 


Figure 20: The red cells denote the aluminum layer, the blue cells denote the air. 


The results are presented in Figure 21. It makes apparent how the temperature of air 
increases between the entry and the exit points. More importantly, it influences the temperature 
distribution of a layer as well. It makes the heat transfer less intensive at the end, thus higher 
temperature difference can occur together with less homogeneous distribution. Consequently, 
the highest temperature can be found at the uttermost point from the entry of the air. 
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Figure 21: The air and the layer temperature distribution in a streamline. The direction of x 
axis is visible in Figure 15. Average airspeed is set to 10 ™, and the heat transfer coefficient is 


62 


W 

m 2 -K' 


3.2.3 Maximum temperature 

One can see the maximum temperature as a function of airspeed in the Figure 22. The 
maximum temperature is lower than the allowed 40 °C if the airspeed is higher than 5™. How¬ 
ever, it is safer to ensure 15™ airspeed, but it is not worth to increase it. The optimum interval 
could be between 10 to 15 m/s. In the further calculations, 10 m/s is used. Unfortunately, that 
speed could be too noisy and too high for clinical applications and may cause other problems 
as well. 

One significant uncertainty is the determination of the heat transfer coefficient. The heat 
convection between the wall and airflow is a very complicated process because of the irregularity 
and uncertainty in the shape of the surface. Hence it is important to investigate the effect of 
changing the heat transfer coefficient. One can see the maximum temperature as a function of 
heat transfer coefficient in the Figure 23. The domain of interest starts at 40 to 80 —Svr. 
Within this domain, the maximum temperature is safely lower than the allowed 40 °C. 

Previously, I mentioned the uncertainty due to the material change that should be investi- 
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Figure 22: Maximum temperature as a function of airspeed. The room temperature is 25 °C. 



Figure 23: The maximum temperature as a function of the heat transfer coefficient. The room 
temperature is set to 25 °C. 
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gated. Now, I handle the thermal conductivity of the absorber layer as a parameter, see Fig. 24 
for details. The domain of interest is 100 77 ^ to 237i which covers all the possibilities for 
aluminium alloys. The maximum temperature variation is found to be 0.3 °C in this interval, 
hence the thermal conductivity of the absorber does not have a big impact on the maximum. 
One of the most likely aluminum alloy for the material of the detector is A11050, which has the 
thermal conductivity of 222^-. The difference between the maximum temperature in case of 
pure aluminum and A11050 is 0.02 °C, which is negligible. 



Figure 24: The maximum temperature as a function of thermal conductivity. 

3.3 Comparison 

The maximum temperature difference in a layer is very important, since the cluster size 
and the noise rate of the ALPIDE are temperature dependent. 1 compared the maximum 
temperature differences for concepts A and B as a function of the airspeed, thermal conductivity 
and heat transfer coefficient. 

One can see the maximum temperature difference as a function of airspeed in case of air 
and water cooling in the Figure 25. As one can see, there is no significant difference in this 
value. If the airspeed is higher than 8 ™ the temperature difference is slightly lower for concept 
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B. Otherwise, the temperature difference is higher than for water cooling. The temperature 
difference in case of concept B, considering 10™ airseed, is 1.6 °C. The maximum temperature 
difference in case of concept A is given as a benchmark, and it is constantly 1.8 °C. In both 
cases, the temperature difference is lower than the allowed 5 °C. 



Figure 25: The temperature difference in a layer in the case of air and water cooling. The 
temperature difference in case of concept A is given as benchmark as it is independent of the 
airspeed. 

The maximum temperature difference depends on the thermal conductivity of the layer, as 
one can see in Figure 26, using 10™ airspeed. One can see that the temperature difference is 
increasing in both cases with decreasing the thermal conductivity. The concept A, using water 
cooling, seems to be more sensitive for this parameter. Considering again the A11050 alloy, the 
temperature difference in the layer is 1.9 °C for concept A and 1.7 °C in case of concept B. 

Changes in the heat transfer coefficient between the layer and the air also has an effect 
in the maximum temperature difference. One can see the outcome in the Figure 27. If the 
heat transfer coefficient is increasing the temperature difference is also increasing. The reason 
behind this effect is the changes in the temperature distribution in the streamline. In Figure 28, 
one can see the temperature distribution using unrealistically low heat transfer coefficient, and 
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Figure 26: Temperature difference in layer in case of air- and watercooling in function of the 
thermal conductivity of the layer. The average airspeed is set to 10 ™ in case of concept B. 

in Figure 29, with unrealistically high heat transfer coefficient. The realistic zone of the heat 
transfer coefficient is between 40 and 80 . In this domain, the maximum temperature 

difference starts from 1.2 °C to 1.8 °C in case of concept B. In case of concept A, the maximum 
temperature difference is not affected by the heat transfer coefficient, so it is only given as a 
benchmark, it is constantly 1.8 °C. 
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Figure 27: The temperature difference in a layer in case of air and water cooling as a function 
of the heat transfer coefficient. The temperature difference in case of concept A is given as a 
benchmark. 



x (m) 


Figure 28: The layer and the air temperature with unrealistically low (15 ) heat transfer 

coefficient, in stream direction. The direction of x axis is visible in Figure 15. The average air 
speed is set to 10 y. 
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Figure 29: The layer and the air temperature with unrealistically high (200 ^ heat transfer 
coefficient, in stream direction. The direction of x axis is visible in Figure 15. The average air 
speed is set to 10 ™. 


36 







4 Cost estimation and environmental effects 

As a result of the previous calculations, I considered that both concept A and concept B 
meets with the requirements of the hadron-tracking calorimeter, so the decision between these 
concepts can be made based on other aspects. One of these aspects is the cost of the cooling 
system, so 1 compare the possible cost of concept A and concept B in this chapter. This 
comparison is a approximation of the costs, so it contains only the cost of the main elements of 
the cooling system. The prices in Table 12 and Table 13 are given as an example for the price 
of the parts. 

4.1 Concept A 

In case of concept "A" I use a water cooling circle. To protect the hadron-tracking calorimeter 
from a possible leak, I need to use water cooling circle, which operates under room pressure. 
If the cooling system operates under room pressure, in case of a leak the air of the room move 
into the tube, not the water move out. To operate under room pressure I need to use a vacuum 
pump. In case of a leak the vacuum pump has to be able to remove all water from the cooling 
circle, so I need the insert a vacuum tank between the cooling circle and the vacuum pump. In 
normal working this tank mainly filled with air. If there is a leak in the cooling circle this tank 
will be filled with the water of the cooling circle. The cooling circle contains a cooler, to cool 
the water under the temperature of the room. The cooler contains a pump, which circulates 
the water in the cooling circle, so I did not calculate the price of the pump separately. The 
cooling circle also contains cooling plates to cool the detector. I only find smaller cooling plates 
then the top or the bottom of the hadron-tracking calorimeter, so I calculated with two cooling 
plate in the top and two cooling plate in the bottom of the detector. The cooling circle also 
contains tubes and fittings. I just estimated price of them, because without an exact plan of 
the cooling system I could not find price for them. The elements of concept A is visible in 
Figure 30. One can see the estimated cost of this solution in Tabic 12. 
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Figure 30: The water cooling circuit of concept A. 


Component 

Main operation data 

Number of components 

Price 

Water cooler [22] 

2400 W, 16.3 4 5 , 2.1 bar 

1 

4090 € 

Vacuum pump [23] 

100 rnbar 

1 

640 € 

Cooling plate [24] 

304.8 mm • 88.9 mm 

4 

110 € 

Vacuum tank [25] 

25 liter 

1 

50 € 

Tubes and fittings 

- 

- 

300 € 

Overall 

- 

- 

5520 € 


Table 12: Cost of concept "A". 


4.2 Concept B 

1 need to calculate the volume flow through the detector and the pressure drop in the 
detector to find a corresponding fan to concept B. 

Pressure drop in case of concept B, if the temperature of the room is T room = 25 °C. The 
equivalent diameter is: 

d ^ = 4 ' 2.L L 2 2 + 2 3 L 3 =L99mm - (14) 

where L 2 = 300 mm and L 3 = 1 mm. The airspeed is w a j r = 10 and the kinematic viscosity 
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(15) 


- 2 

is v = 1.58 • 1CT 5 The Reynolds-number of the airflow is: 

Re = Wair ' dcq = 1250, 
v 

which means the airflow is laminar between the detector layers. The pressure drop can be 
calculated with this equation 


Ap 


3 • Lin • pw air 



438 Pa 


(16) 


based on Newton’s law of viscosity, where p — 1.19 and Li = 200 mm. The volume flow 
between the layers is 


Q ~ n layers • L3 • L2 • W a j r 


0.105 


m 



(17) 


where n layers = 35. 

1 need to have a fun in the cooling system to generate airflow in the detector. The fan is as 
big as the hadron-tracking calorimeter and it is possible the fan generate vibrations, so it has 
to operate separately from the detector connected with a tube. 1 also need an air filter before 
the detector, to protect the hadron-tracking calorimeter from the dust in the air. The elements 
of the cooling system is visible in Figure 31. 



Figure 31: The elements of B cooling system concept. 
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Component 

Main operation data 

Number of components 

Price 

Fan [26] 

950 Pa, 400 ^,485W 

1 

890 € 

Tube [27] 

1=10 m, d=250 mm 

1 

50 € 

Fittings and filter 

- 

- 

300 € 

Overall 

- 

- 

1240 € 


Table 13: Cost of concept "B". 


4.3 Comparison 

As a result of the cost estimation one can conclude that the concept B is cheaper than 
concept A. As environmental effect one can consider that concept B use less energy to operate 
than concept A, because the water cooler needs more energy to cool the water under room 
temperature, than the fan needs to transfer air through the detector. 

It is possible to use a concept like A, but use thermoelectric cooler and heat sinks with 
heat pipe and funs to cool the top and the bottom of the detector. This solution results same 
temperature distribution than concept A, but it can be cheaper, smaller than concept A. In 
case of this concept it is also possible to cool the detector under room temperature and the 
maximum temperature also adjustable. This concept can be more advantageous than concept 
A in case of beam test, but concept A is more suitable for clinical use, because in this case the 
generated heat is transported out of the treatment room. 
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5 Finite element simulation of concept A 


My next work in the research team was the more accurate calculation of the temperature 
distribution in case of concept A. I decided to check the effect of contact thermal resistances and 
inhomogeneous load on the temperature distribution. I took into account the heat generation of 
the ALPIDEs as a heat flux in the contact surfaces of them. I calculated with a two dimensional 
model (Figure 32) when I calculate the effect of the contact resistances, and I calculate with a 
three dimensional model (Figure 11) when I calculated the effect of the inhomogeneous load. I 
used Aauo50 = 222 as the thermal conductivity of the absorber layers in both calculation. 

5.1 The effect of contact thermal resistances 

Firstly I calculated with 4 mm thick absorber layer. In reality the layer will be build from 2 
mm thick absorber layer in the middle, and two 1 mm thick layer in both side. The ALPIDEs 
will be glued into the 1 mm thick plates, as one can see in 32. There will be contact resistances 
between this plates and between the cooling plate and this plates. I calculate the contact 
resistance between the plates as the resistance of the average air gap between the layers, based 
on [28]. I calculate the average air gap between two plate with = R z equation. The 

aluminum plate is cold rolled aluminum, so the mean roughness depth is R z = 12.5 fin i. The 
average distance is <5i = 12.5 /jrn. Between the cooling plate and the edge of the absolber 
plates I calculated with S 2 = 200 fix n, because of the uncertainty of the positioning of the 
plates. It is necessary to use thermal conductive grease between the cooling plate and the edge 
of the absolber plates. I considered A a i r = 0.0261 ^ 7 ^ as the thermal conductivity of air and 
^conductive = 8.5 [29] as the thermal conductivity of a typical thermal conductive paste. 

The thermal resistance between two layer is Rq = = 4.79 • 10 4 ni ^ 4v and between the 

cooling plate and the edge of the absorber layers is R 2 = t—— = 2.35 • 10 5 m ,,f . 

^conductive 
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Figure 32: The model of the layer with thermal resistances 


One can see the model of the layer and the contact resistances in Figure 32 and the result 
of the two dimensional simulation in Figure 33 and Figure 34. 1 calculated the maximum 
temperature difference in that part of the surface of the absorber layer, where ALPIDEs are 
glued. 
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Figure 33: The effect of the thermal resistance between layers in the maximum temperature 
difference. 



Thermal insulance 


Figure 34: The effect of the thermal resistance between the cooling plate and absorber the 
layers in the maximum temperature difference. 


I have checked the mesh dependence of the solution in case when Ri is 4.79 • 10 4 and 

_ 2 I/ 

R.2 is 2.35 • 10 0 ni ^r\ I checked the maximum temperature difference (AT max ) in y=2 mm 
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line, which is visible in Table 14. I found that the mesh with 0.5 mm edge length seems to be 
the most accurate. The maximum temperature difference between 1 mm edge length and 0.5 
mm edge length is small compared with the measured temperature differences, so the solution 
is mesh independent. 


L in ax in x direction (mm) 

L in ax in y direction (mm) 

Number of elements 

ATmax (°C) 

1 

1 

800 

- 

0.5 

0.5 

3200 

0.017 

2 

0.5 

824 

- 

1 

0.25 

3216 

0.104 


Table 14: Outcome of the mesh dependence test, where AT max is the maximum difference of 
temperature compared with the previous mesh. L max is the maximum edge length in the given 
direction. 

I considered based on Figure 33 and Figure 34 that the contact thermal resistance between 
the layers has minimal effect on the temperature difference, and the contact thermal resis¬ 
tance between the cooling plate and the layers has only marginal effect on the temperature 
distribution. 

5.2 The effect of inhomogeneous load 

I obtained the temperature distribution in the case, if only the middle of the detector is 
used. I expect the maximum temperature difference will be smaller than in case of full load, 
but it is interesting how much it will be, because the partial load will be a typical usage, if we 
have to take image from a smaller part of the human body then the size of the full detector 
area. 

I used P s t a ndby = 60 mW heat generation in every ALPIDE detector, which was inactive in 
the measurement, which is the power consumption of waiting for hits. I used P ac tive = 185 mW 
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heat generation in every active ALPIDE, which is the heat generation working with maximum 
readout capacity. 

I obtained the temperature distribution in an x direction path, defined by z=0 mm and 
y=0.002 mm. The temperature distribution in this path is visible compared with the temper¬ 
ature distribution of full load in Figure 35. 



x [m] 


Figure 35: The temperature distribution in case of full and partial load in path z=0 mm and 
y=0.002 mm. 

One can see the temperature in a path, defined by x=0 mm and y=0.002 mm, in Figure 36. 
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Figure 36: The temperature distribution in case of full and partial load in path x=0 mm and 
y=0.002 mm. 


The result of the mesh dependence test is visible in Table 15. 


Maximum edge size (mm) 

Number of elements 

AT max (°C) full load 

AT max (°C) full load 

7.5 

1234 

- 

- 

5 

2400 

0.171 

0.077 

2.5 

19200 

0.086 

0.038 


Table 15: Outcome of the mesh dependence test in case of full and partial load, where the 
AT max is the maximum temperature difference compared with the previous mesh. 

The temperature difference was 2.1 °C in case of full load and 1.1 °C in case of partial load. 
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6 Calculation of transient behavior 


In this subsection I focus on how much time does the heat up takes, because it effects both 
the clinical use, both the measurements during the development of the detector. I used the 
same geometry than I used in Chapter 3.1, so I used a one dimensional model. I calculated 
with the thermal conductivity of A11050, which is 222 The heat transfer is described by 
the following equation: 


dT _ <9 2 T q v 

dt a 'dx 2+ p- c 


(18) 


where a = 

Initial condition: 


T(x, 0) = T 0 = 25 °C. 


(19) 


Boundary conditions: 


T(0,t) = Ti = 25 °C (20) 


and 


T(0,t) = T 2 = 25 °C. (21) 


The temperature can be calculated as the sum of steady state and the transient temperature: 


T(x,t) = T s (x) + T h (x,t). 


( 22 ) 


The steady state temperature is [20]: 


T 8 (x)=-^~+C 1 -x + C 2 , 


where C 2 = Ti and Ci = T2 L Tl + x ' ?• The transient temperature is [20]: 


(23) 


Tp(x,t) = £ Dn • e 


,-A n -t 


Sill 


n=l 


n ■ 7T 


L 


x 


(24) 


1 a fn-Tr \ 2 1 ~rx fn J (T( x > 0 )-T s (x))-sin(Sp E -x)dx 

where A n = a • W and D n = Jo v v V —^ 

V L / / 0 L sin 2 (a^-x)dx 

One can see the temperature of the middle of the detector in function of time in Figure 37, 


and the parameters, which was used to the calculation, in Table 16. 
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Figure 37: The temperature of the middle point of the detector in function of time. 


Parameter 

Value 

Initial temperature of the layer (Ti) 

25 °C 

Temperature of the edge of the layer (Ti and T 2 ) 

25 °C 

Thermal conductivity of the layer (A) 

222 \ 
rii-K 

Mass density of the layer (p) 

2710 H 

m° 

Specific heat of the layer (c) 

434 yh 

Equivalent volumetric heat generation (q v ) 

8.3 • 10 4 ^ 

m° 

Maximal item of Fourier series (n max ) 

50 


Table 16: Parameters used in the heat up calculation. 

It takes 134 s, while the temperature of the middle point of the layer reach ±0.1 °C zone 
around the steady state temperature of the middle point. 
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7 Looking Ahead 


In this work I analysed the temperature distribution of a layer of the hadron-tracking 
calorimeter. The detector is going to contain two layer before the calorimeter, which mea¬ 
sure the path of the particles. This layers are going to contains a carbon composite layer in 
the middle instead of an aluminum absorber layer. This carbon composite layer is going to 
work as a support structure, and it is thinner than aluminum absorber layer. The thermal 
conductivity of carbon composite is anisotropic, it is reasonable in longitudinal direction, but 
low in perpendicular direction. It is essential to analyze the temperature distribution in this 
layers, and design new cooling concepts for these layers, if it is necessary. 

Another important engineering task in the development of proton computer tomograph 
is the comparison of the temperature distribution between layers. This task requires more 
information about the power consumption of the tracking detectors and the number of particles, 
which goes through them. It is possible to reach this information with test beam measurements, 
so one of the next engineering task should be the measurement of the electric power consumption 
in case of different particle fluxes. 
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8 Summary 


Cancer started to become the leading cause of death in the developed world, as it is re¬ 
sponsible for 25% of all death in Hungary. Mainly, there are three ways of treatment: surgery, 
radiotherapy, and chemotherapy. Hadron therapy is a novel treatment in radiotherapy. It is 
beneficial because it generates less radiation in healthy tissues than X-ray radiotherapy, thus 
it leads to fewer side effects, and it allows a higher daily dose, which increases the effectiveness 
of treatment. 

Hadron therapy requires an accurate three-dimensional map about the electron density of 
the body. The most accurate solution to obtain the map uses hadrons, however, it requires the 
development of a proton computer tomograph (pCT). Our group, the Bergen pCT Collabora¬ 
tion, aims to develop a pCT that can meet the requirements of the clinical use. The detector 
of the pCT is built of silicon pixel tracking detectors and aluminum absorber layers. 

In the collaboration, my role was to investigate the temperature distribution of the detector, 
which is essential because of the accuracy of the energy measurement based on the temperature 
of the detector. It is indispensable in the development of the pCT detector. 

I applied the lumped capacitance method to estimate the thermal behaviour of the detector, 
containing 35 absorber layers. This simple analysis shows the necessity for cooling. Hence, two 
cooling concepts are elaborated in this work. In concept A, water cools the edges of the 
aluminum absorber layers. In concept B, the air is circulated between the detector layers. My 
goal was to quantify the differences in the steady-state temperature distribution and compare 
the concepts to each other. 

Firstly I calculated with analytic methods the steady state temperature distribution in 
case of both concepts. Both of the them satisfies the first requirement, which is to keep the 
maximum temperature lower than 40 °C. The allowed temperature difference is 5 °C. Both of 
the concepts meets this requirement also, but concept B has a bit better performance as the 
temperature difference is 1.6 °C in this case, compared with 1.8 °C in case of concept A. 

Secondly I estimated the cost of both concept. I considered concept A is more expensive 
than concept B. It is possible to use a third solution, which cools the top and the bottom of the 
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hadron-tracking calorimeter, like concept A, but contains thermoelectric coolers, heat pipes, 
heat sinks and funs instead of a water cooling circuit. This concept can be smaller and cheaper 
than concept A, but it has the advantageous properties of concept A. 

Thirdly I used finite element simulations to determinate the effect of the contact thermal 
resistances and partial load in the temperature distribution, in case of concept A. I obtained 
that a reasonable contact thermal resistance has small effect on the temperature distribution, 
as it cause maximum 0.2 °C increase in the temperature difference. In case of partial load the 
temperature difference is decrease, as it was expected. I calculated the time, which is required 
to reach steady state temperature distribution. It takes 134 s, which is suitable for clinical use 
and test beam measurements. 
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9 Osszefoglalo 

Napjainkban a rak az egyik vezeto halalokka kezd vallni a fejledt vilagban, Magyarorszagon 
a tcljes clhalalozas 25%-aert felelos. A raknak harom fele fo gyogymodja van: az elso a sebeszet, 
a masodik a sugarterapia es a harmadik a kemoterapia. A hadronterapia egy iij sugarterapias 
kezeles. Eldnyds, mert kevesebb ionizaciot okoz az egeszseges szovetekben mint a hagyomanyos, 
rontgen sugarzassal vegzet sugarterapia. Ennek kdszdnhetden kevesebb mellekhatassal jar, es 
nagyobb napi dozist tesz lchetdve, ezzel novelve a kezeles hatasossagat. 

A hadronterapiahoz sziikseges egy pontos terbcli kepet kesziteni a paciens testenek elektron- 
surusegerol. A legpontosabb megoldas, ha ezt a kepet hadronok hasznalataval keszitjiik el, de 
ez egy proton computer tomograf (pCT) kifejleszteset igenyli. A Bergen pCT egyuttmukodes 
cclja, hogy kifejlesszen egy ilyen eszkozt, ami a klinikai felhasznalasi szempontoknak eleget tesz. 
A pCT detektora sziliciumpixel-nyomkoveto detektorokbol es aluminum abszorber retegekbol 
epiil fel. 

Az egyiittmukodesben betoltott szerepem a detektor homerseklet-eloszlasanak meghataro- 
zasara iranyult. Ennek vizsgalata clcngedhetetlen a pontos meresek elvegzesehez. Elbszor, 
egy egyszeru, koncentralt parameteru modellcl becsiiltem a 35 retegbol alio detektor termikus 
viselkedeset. Ebbol egyertelmuen latni, hogy a hutes elengedhetetlen. Ket fele hutesi koncepciot 
dolgoztam ki. Az “A” koncepcio eseten vizzel hutjom az abszorber retegek szelcit, mig a “B” 
koncepcio eseten levegot keringtetek a retegek kozott. Celom a kialakulo, allandosult allapotu 
homerseklet-eloszlasok meghatarozasa es a ket koncepcio osszehasonlitasa. 

Eloszor analitikusan meghataroztam az allandosult homersekleteloszlast mindket koncepcio 
eseteben. A megengedett maximum homerseklet a detekteorban 40 °C, amit mindket koncep¬ 
cio teljesit. A megengedett hdmerseklet-kulonbseg 5 °C, amit ugyszinten mindket koncepcio 
teljesit, de a B koncepcio ilyen szempontbol clonyosebb, mivel ebben az esetben a homerseklet- 
kulonbseg csak 1.6 °C, amig az A koncepcio eseteben 1.8 °C. 

Masodszor megbecsultem a ket koncepcio lehetseges koltsegeit. Eredmenyul azt kaptam, 
hogy az A koncepcio lenyegesen dragabb. Lehetseges egy harmadik koncepcio alkalmazas is, 
amcly eseteben ugyanugy a detektor also es felso fclulctet hiitjuk, mint az A koncepcio es- 
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eteben, azonban vfzhutes helyett Peltier-elemek, hocsovek, bordak es ventillatorok segftsegevcl 
valositjuk meg a hutest. Ez a koncepcio kisebb es olcsobb megoldast jclcnthet az A koncep- 
cional, azonban rendelkezik annak az elonyos tulajdonsagaival. 

Harmadszor vegeselemes-szimulaciok segftsegevcl vizsgaltam a kontakt-hdellcnallasok es az 
reszleges terheles hatasat a homerseklet eloszlasra az A koncepcio eseteben. Eredmenyiil azt 
kaptam, hogy a realis kontakt-hdellcnallasok maximum 0.2 °C hbmersekletkulonbseg-novekedest 
okoz. Reszleges terheles eseten a maximalis hdmerseklet-kiilonbseg kisebb, mint teljes ter¬ 
heles eseten, ahogyan ezt elore vartuk. Analitikus szamitassal meghataroztam az allando- 
sult homersekleteloszlas kialakulasahoz sziikseges idot. Ez 134 masodpercnek adodott, ami 
megfelcld, mind az orvosi fclhasznalas, mind a tesztmeresek szamara. 
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